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Abstract We present the results of analytical modelling of fast-mode magne- 
tohydrodynamic wave propagation near a 2D magnetic null point. We consider 
both a linear wave and a weak shock and analyse their behaviour in cold and 
warm plasmas. We apply the nonlinear geometrical acoustics method based on 
the Wentzel-Kramers-Brillouin approximation. We calculate the wave ampli- 
tude, using the ray approximation and the laws of solitary shock wave damping. 
We find that a complex caustic is formed around the null point. Plasma heating 
is distributed in space and occurs at a caustic as well as near the null point 
due to substantial nonlinear damping of the shock wave. The shock wave passes 
through the null point even in a cold plasma. The complex shape of the wave 
front can be explained by the caustic pattern. 

Keywords: Magnctohydrodynamics; Waves, Propagation; Waves, Shock; Heat- 
ing, Coronal 

1. Introduction 

Various eruptive processes in the solar atmosphere as well as convective mo- 
tions beneath the photosphere generate magnetohydrodynamic (MHD) waves. 
Propagation of waves in an inhomogencous plasma can result in wave energy 
dissipation, which is known to contribute into the general energy balance of the 
solar corona. Behaviour of MHD waves is governed by the topology of mag- 
netic field and gradients of MHD quantities. Singularities of the magnetic field 
topology are null points, at which the magnetic field magnitude is equal to zero. 

Dissipation of the wave energy is believed to occur most efficiently near mag- 
netic null points. In the solar corona, where the plasma beta is much less than 
unity, a fast-mode MHD wave is subject to significant deceleration at nulls, which 
therefore are the regions of the accumulation of the wave energy. The amplitude 
of the wave and gradients of all MHD quantities increase considerably due to 
the Alfven speed decrease as well as the wave front convergence in the vicinity 
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of the null point, where the wave becomes nonlinear. The wave profile steepens 
and the initially linear wave transforms into a shock. Another important effect 
of the increase of the gradients is the appearance of steep spikes of the electric 
current density in the vicinity of the null point, as it was demonstrated by 
IFoullon et aLI(|2005p and modelled numerically by |McLaughlin and Hood| (|2"0"0"3[) 
and iNakariak ov et aTl (|2006|) . These circumstances result in rapid conversion of 
the wave energy into heat. 

Topological properties of magnetic null points have been studied well (see, 
e.g., the review by |Longcop"e[ I2005P and the behaviour of waves in the vicin- 
ity of a null point has been investigated by many authors (see the review by 



McLaughlin, Hood, and De Moortel| EMTjl . |McLaughlin and Hood] (j2"0"M)) con- 
sider the incidence of a linear fast-mode wave on a 2D null point both numerically 
and analytically using the Wentzel-Kramers-Brillouin (WKB) approximation. 
The authors neglect the plasma pressure against magnetic pressure and as- 
sign the sound speed to be zero. Their analysis shows that ray paths along 
which the wave travels get curved towards the null point due to refraction. 
The wave envelops the null point and its energy is converted to heat. Then 
McLaug hlin and Hood| (|2006ap consider the wave incidence on a 2D null point 
formed by two magnetic dipoles. In this case, only a portion of the wave and its 
energy are captured by the null point. The rest propagates away. 

The assumption of zero sound speed is a significant restriction for the problem 
of interest. The solar coronal plasma has a temperature of about 1.5 x 10 6 K, 
which gives a sound speed of about 180 km s _1 . There is an area around the 
null point, in which the sound speed c is greater than the Alfven speed Va- This 
considerably affects wave propagation. The wave is able to pass through the null 
point in this case. A numerical modelling and analytical analysis of this problem 
is carried out, in the linear approximation, by |McLaughlin and Hood] (|2006b[l . 
In their numerical study, they also find conversion of the wave modes when the 
wave crosses the Va = c layer. 

Further, [McLaugh lin, Ferguson, and Hood|(|2008p consider the wave propaga- 
tion near a 3D magnetic null point under the assumption of cold plasma, using 
the WKB approach. This investigation confirms capturing the wave by the null 
point as well as plasma heating at that place. We note that in the analytical 
studies mentioned, the authors model only the wave front propagation and do 
not examine the wave amplitude evolution. 

The propagation of a nonlinear wave near a magnetic null point is inves- 
tigated by |McLaughlin et al . (2009) for the magnetic reconnection problem. 



Their numerical modelling shows a complex pattern including the formation 
of a current sheet as well as the generation of secondary shock waves. Using 
a similar numerical approach, IGruszecki et al.\ (|201ip study the nonlinear wave 
propagation to understand the triggering of the sympathetic solar flares. 

The aim of this paper is to analyse analytically nonlinear effects of the weak 
fast-mode MHD shock wave propagation near a magnetic null point and calculate 
the wave amplitude evolution for the coronal heating problem. We do not study 
mode conversion as well as the transformation of the null point into a current 
sheet by the shock wave, since the ray method is not able to take this into 
account. The paper has the following outline. We consider the linear wave case 
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in Section [3J then the case of a shock wave in Section [3] and discuss the results 
obtained. Section 2] contains concluding remarks. In the Appendices we give 
bulky mathematical details of our investigation. 



2. Linear wave 



In this study we explore a two-dimensional magnetic null point, which is lo- 
cated at the point of origin and specified in Cartesian coordinates (x, y, z) as 
B = (x, 0,— z) (see, e.g.. IParnell et a/lll996p . In this model, the magnetic field 
magnitude grows linearly with the distance from the null point. The plasma 
density p and temperature T are assumed to be constant, so the sound speed 
c = J 2~fR gas T/ ffiMfj is also constant and the Alfven speed Va = B/yJAirp 
distribution is governed only by that of the magnetic field. Here Mh is the 
molar mass of hydrogen, m the average atomic weight of an ion, 7 the adiabatic 
index, and R gas the gas constant. The Alfven speed vanishes completely at the 
null point and grows to infinity away from it. 



2.1. Cold plasma 



At first we consider the case of a cold plasma, with the sound speed c being zero 
overall. Let a linear plane fast-mode MHD wave propagate towards the magnetic 
null point. To calculate its propagation, we apply the geometrical acoustics 
method (WKB approach). In this method, a solution is found in the form of 
A (r, t) e 4 *( r,t ) where A(r,t) is the wave amplitude, and ^ (r, t) is the eikonal, 
depending on both coordinates r = (x, z) and time t. By substituting this rep- 
resentation for wave perturbations into the set of linearised equations of ideal 
magnctohydrodynamics, one can obtain a Hamilton-Jacobi partial differential 
equation for the eikonal of fast-mode MHD waves (see, e.o,. IUralova and Uralovl 
IT3MD : 

a* 

— + a|grad*| = 0, 

where a = 1/2 {^/c 2 + V\+ 2cV A cos a + y 1 c? + V%- 2cV A cos aj is the fast- 
mode velocity in plasma, a the angle between the direction of the wave prop- 
agation and the magnetic field vector. Solving the equation with the method 
of characteristics gives the set of ray equations, which are a set of ordinary 
differential equations. In Appcndix[X]wc provide this set of equations (see Equa- 
tions Q). |McLaughlin and Hoo d (2004) succeeded in finding an exact analytical 
solution of the set for the magnetic field model under consideration: 

_itt t / X , : . Xq \ _£Q t / Xq . Xq 

X = e r o Xq COS 1 + Zq Sin 1 , Z — e ' Zq COS 1 — Xq Sin 1 

\ Tq Tq ) \ Tq Tq 

(1) 

where ro = y x o + 2 0' ano - and zq are the initial values of a ray path. 

The null point attracts the wave, so all the ray paths tend to the origin where 
the null point is located. Figure [1] shows the ray pattern. In both the figures and 
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Figure 1. The propagation of a linear wave near a null point in a cold plasma. The thin lines 
represent the rays and the thick one marks the wave front at moment t = 1. At the initial 
moment, the plane wave front was located at z = 1 and — 1 < x < 1. The axes are shown in 
the normalised units. 



calculations throughout the paper, we use normalised units. The normalising 
values are the characteristic spatial scale of the problem, and the magnetic field 
strength and the Alfven speed at the distance of the characteristic spatial scale 
from the null point. The normalising time value is then set as the ratio of the 
characteristic spatial scale to the characteristic Alfven speed. 

However, it is essential not only to find the wave geometry, but also to 
calculate the wave intensity. Geometrical acoustics allows the wave amplitude 
variation to be calculated. 

In linear geometrical acoustics, the energy flux of a disturbance travelling at 
group velocity q in a stationary medium is directed along the rays, and its mag- 
nitude is conserved within a ray tube |Blokhmtscv, 1981), div(Aeq) = 0, where 
As = p(u 2 + v 2 ) is the average density of the disturbance energy, p the undis- 
turbed plasma density, and u, v the plasma velocity components along the normal 
k to the wave front and across it, respectively. Taking into account the relation 
between the plasma velocity components ( |Kulikovsky and Lyubimov, 2005[ ), p, = 
vju = (l — c 2 / a 2 ) cot a, it is possible to relate the variation of the wave ampli- 
tude to the cross section dS of the ray tube formed by a bundle of rays: 



dSqpu 2 (l + /i 2 ) = const. 



(2) 



We calculate the cross section, using the Jacobian of the transformation from 



Cartesian coordinates to ray ones ( Kravtsov and Orlov, 1990 ). The volume el- 
ement dW of the ray tube is expressed in terms of the ray coordinates (77, t) 



as 



dW = dxdz = D (t) drjdt, D (t) 



dx/dn dx/dt 
dz/dn dz/dt 



where D(t) is the Jacobian. Then for the cross section of the ray tube, dS, we 
have 



dS=~-- = ^dn 



dW 
da 
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Figure 2. Evolution of the linear wave amplitude along the rays with xrj = 0.5 (a) and xq = 1 
(b) (marked as the dashed lines in Figure[T]| in the cold plasma case. The dashed lines represent 
the plasma velocity component u along the normal to the wave front and the solid ones show 
the magnitude of the plasma velocity. The horizontal axes are shown in the normalised time 
units, the vertical axes are shown in the units of the initial amplitude. 



where a is the ray length. Let the ray coordinate r\ be the initial value Xq. Using 
Equations (JTJ) and @ we obtain the expressions for the wave amplitude A along 
the ray with the initial values Xq and zq (Aq is the initial amplitude): 

J-P— e&, (3) 
\ Z t + r W 

if we assume the wave amplitude A to be the component u of the plasma velocity 
vector along the normal to the wave front, or 

A(t)/A = ,r^ e- 4 , (4) 
V z t + r 

if we assume A to be the magnitude of the plasma velocity vector, A = y/u 2 + v 2 . 

Figure [5] shows time evolution of the wave amplitude along two different rays. 
The amplitude plots in the paper are shown in the units of the initial amplitude. 
The plasma velocity component u oscillates and drops to zero from time to time, 
while the magnitude of the plasma velocity grows exponentially. Such behaviour 
of the component u is due to changing the direction of the wave propagation. 
The wave front element travels around the null point and sometimes the wave 
front normal k becomes collinear to the magnetic field B, i.e. a — or tt. At 
these points, the component u disappears but the component v becomes equal 
to the magnitude of the plasma velocity vector at that moment. There is also 
the opposite case: when the a = ±7r/2 the component v drops to zero. So the 
solid line marking the magnitude of the plasma velocity is an envelope line. 
Note that in Figure [5] we plot a relative variation of the amplitude and adjust 
the magnitude of the plasma velocity vector to be unity at the initial moment. It 
means that for the component u we actually use Equation ([3]) divided by ?"o/|xo|. 

2.2. Warm plasma 

Now we consider the case of a warm plasma with non-zero sound speed. Let 
points where the Alfven speed is equal to the sound speed be named the Va = c 



A(t)/A 



x cos2^t + z sin 2^4 

u r u ra_ 

X 
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layer (see, e.g., |McLaughlin and Hood[ I2006b|) . In the problem of interest, the 
case of a warm plasma corresponds to the situation when the wave propagates 
near the Va = c layer, while in the case of a cold plasma the wave propagation 
is considered far from it. The ray equations become more complex and we have 
to integrate them numerically. To calculate the Jacobian and find the wave 
amplitude, we use a method based on numerical integration of the so-called 
adjoint set of equations. Solving the adjoint set allows the quantities dx/drj and 
dz/drj to be calculated. For more details, we refer the reader to the paper by 
| Afanasyev and Uralov|P01ip . Appendix lAl contains the set of ray equations (see 
Equations (O) as well as the adjoint set (see Equations ([5])), which we solve to 
model the wave propagation. 

Like McLaughlin and Hood (|2006b|) we find that the wave passes through 



the null point. Figure [3] shows the wave propagation in this case. We study 
behaviour of the wave front, which was located at z = 1 and —2 < x < 2 at 
the initial moment, however we plot the rays outgoing from a central part of the 
wave front for clearing the ray pattern. Unlike the cold plasma case, ray paths 
begin to intersect near the top boundary of the Va = c layer. This results in the 
considerable distortion of the wave front as well as the subsequent formation of 
a caustic. 

By definition, a caustic is an envelope of a family of rays, i.e. rays approach 
each other at points of the caustic ( |Kravtsov and Orlov, 1990] ), and therefore 
a caustic can be identified in the ray pattern as a ray crowding. The Jacobian 
and the ray-tube cross sections drop to zero at the caustic points and the wave 
amplitude grows indefinitely. Figure [4] shows evolution of the wave amplitude 
along two rays outlined by red in Figure [3] Marking points of the ray paths 
where the amplitude tends to infinity for the first time (since this can occur 
twice along some rays during the calculation time in our modelling) we can 
obtain the caustic line. The bottom panel in Figure [3] shows the inner part of 
the caustic, associated with a central part of the wave. The rest of the caustic 
encompasses the null point, bending away from it. The outer parts of the caustic 
form if the initial front is rather extensive. Besides, if the calculation time is also 
rather long it is possible to identify secondary (repeated) caustics. For instance, 
one of them can be seen in the top panel of Figure [3] as a cusp at the lower part 
of the plot. Note that in the linear cold plasma case there are no caustics in the 
neighbourhood of the null point. 

The singular behaviour of the wave amplitude indicates that the ray approx- 
imation does not work at the caustic. Indeed, more accurate considerations show 
that the amplitude grows but it always remains finite ( |Kravtsov and Orlov, 1990| .[ 
Nevertheless, caustics give important information on the amplitude distribution 
along the wave front. At every moment, the wave amplitude increases signifi- 
cantly at a few points (the caustic does not "fire up" simultaneously throughout 
its extent). Caustics in the ray pattern represent the places of the most efficient 
plasma heating associated with the wave. The energy of the wave accumulates 
in the neighbourhood of the null point, but this accumulation is distributed in 
space, unlike the cold plasma case. Our calculations show that the size of the 
Va — c layer may be chosen as the characteristic spatial scale of the heated 
region around the null point. 
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Figure 3. The propagation of a linear wave near a null point in a warm plasma: the ray 
pattern (top), and the wave front and the caustic (bottom). In the top panel, blue denotes the 
location of the Va = c layer and the null point, black and red — the rays. In the bottom panel, 
the solid lines denote the wave front at moment t = 2.4 (thick) and the caustic (thin), and the 
dashed line denotes the Va = c layer. At the initial moment, the plane wave front was located 
at zq = 1 and —2 < x < 2. The axes are shown in the normalised units. 
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Figure 4. Evolution of the linear wave amplitude along the rays with xo = 0.5 (a) and 
xq = 1 (b) (red lines in Figure [3j in the warm plasma case. The dashed lines represent the 
plasma velocity component u along the normal to the wave front and the solid ones show the 
magnitude of the plasma velocity. The horizontal axes arc shown in the normalised time units, 
the vertical axes are shown in the units of the initial amplitude. 



Analysing the caustic pattern we can also readily explain the complicated 
surface of the wave front mentioned in McLaughlin and Hood (2006b). A break 
of the wave front occurs at the caustic (see Figure[3]), so "small triangular shapes" 
of the wave front arise at the caustic cusps. Vertices of these triangles travel along 
the caustic lines as the wave front propagates. 



3. Shock wave 

As seen in Section 12.11 the wave amplitude grows exponentially as the wave 
approaches the null point. This results in the nonlinear steepening of the wave 
front and therefore a shock wave formation. So we should take into account 
nonlinear properties of shock waves and investigate how they affect the wave 
propagation. 

Firstly, a nonlinear plane disturbance in an ideal homogeneous medium prop- 
agates at speed that exceeds the fast-mode velocity and depends on the dis- 
turbance amplitude. In a nonlinear acoustics approximation, a fast-mode wave 
element with plasma velocity component u normal to the front moves at a speed 
of a + ku, where k = (1/a) (d(pa) /dp) is the numerical coefficient depending 
both on the plasma beta and the angle a between the wave front normal and 
the magnetic field (see Appendix [A}. The fact that each element of the wave 
travels at its own speed causes the wave profile deformation and the appearance 
of a shock wave. In turn, a weak shock moves at a speed of a + KU a h/2, where 
U s h is the jump of the plasma velocity component u in the shock. 

Secondly, the nonlinear wave amplitude undergoes additional damping asso- 
ciated with the energy dissipation in the discontinuity. On the one hand, the 
shock wave amplitude increases due to the convergence of rays as well as the 
Alfven speed decrease. On the other hand, the amplitude decreases due to the 
nonlinear damping. 

The above-mentioned basic properties of shock waves arc taken into account 
in the method of nonlinear geometrical acoustics. It allows one to calculate the 
propagation of disturbances with small but finite amplitudes (i.e. weak shock 
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Figure 5. The propagation of a shock wave near a null point in a cold plasma. Blue colour 
denotes the positions of the wave front at moments t = 2.5 and t = 3.0, green - the null point, 
black and red - the rays. At the initial moment, the plane wave front was located at zq = 1 
and — 1 < x < 1. Panel b shows a zoom on the null point vicinity. The axes are shown in the 
normalised units. 




Figure 6. Evolution of the shock wave amplitude along the ray with xq = 0.5 in the cold 
plasma case. The thick lines correspond to the amplitude of the shock wave and the thin ones 
correspond to the amplitude variation due to only the convergence of rays and the Va decrease 
(without dissipation in the discontinuity). The dashed lines denote the jump U a h of the plasma 
velocity component u in the shock wave, the solid ones - that of the plasma velocity magnitude. 
The horizontal axis is shown in the normalised time units, the vertical axis is shown in the 
units of the initial amplitude. 



waves) through an inhomogeneous medium. To familiarise the readers with de- 
tails of the method used, we refer them to the papers by |Afanasyev and Uralov| 
(|201ip and lUralova and Uralovl (|1994|) . Besides, in Appendix lAl we give the set 
of nonlinear ray equations (see Equations ([5])). which we solve to model the 
weak shock wave propagation, and in Appendix [B] we present derivation of the 
nonlinear ray equations. 

Figure [5] shows the propagation of a weak shock wave in a cold plasma. 
Like the linear cold plasma case, the wave is captured by the null point and 
wraps around it. However, the propagation speed at the null point is non-zero 
owing to nonlinearity and the shock wave passes through the null point (cf. 
|McLaughlin et ql] I2009[) like a linear wave in a plasma with the small sound 
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Figure 7. The propagation of a shock wave near a null point in a warm plasma. The dashed 
line denotes the location of the Va = c layer, the solid ones are the rays. At the initial 
moment, the plane wave front was located at zrj = 1 arid — 1 < x < 1. The axes are shown in 
the normalised units. 

speed. Therefore, in this case a caustic also arises in the small neighbourhood 
of the null point. This can be seen in Figure [SJd, which shows the propagation 
pattern in a small area around the null point. The caustic is seen as the sharp 
border of a ray crowding. The amplitude of the wave tends to infinity at the 
caustic. This indicates again the places of significant plasma heating. 

The wave front elements being close to the null point have higher amplitudes 
and therefore they travel faster. This results in forming a bulge in the wave front, 
which is visible in Figure [SJd underneath the null point. 

Figure [5] shows the shock wave amplitude. It decays considerably due to 
dissipation of the wave energy in the shock front. For comparison we also plot 
variation of the wave amplitude calculated with the assumption that there was 
no dissipation in the shock front. These plots demonstrate that a significant part 
of the shock wave energy is spent on plasma heating. 

We have also considered the propagation of a weak shock wave in a warm 
plasma. In this case both nonlincarity of the wave and the non-zero sound speed 
in a plasma work in the same way. The shock wave propagates through the null 
point and a caustic is formed. The ray pattern is shown in Figure [71 The caustic 
is quite similar to that in Figure [3] for the linear warm plasma case, however 
it is shifted in the direction of the wave incidence. An effect of the non-zero 
sound speed appears to predominate over that of the wave nonlinearity. This is 
not unusual if we keep in mind that the nonlinear geometrical acoustics method 
should be applied only for weak shock waves. We also note that two caustic cusps 
are located at the Va = c layer. 
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4. Discussion and Conclusions 

We have applied the nonlinear geometrical acoustics method to model the prop- 
agation of a fast- mode MHD wave near a 2D magnetic null point. We have 
considered both a linear wave and a weak shock and calculated their behaviour 
in cold and warm plasmas. We have also calculated the wave amplitude, using 
the ray approximation and the laws of solitary shock wave damping. 

The nonlinear ray method allows one to trace the wave front and calculate the 
jumps of MHD quantities at the front only, but it does not give information on 
the downstream plasma flow evolution, since irreversible processes occur behind 
the shock front. In the propagation model used taking into account the refraction 
effect (i.e. the wrapping effect), the wave front can encounter plasma "upstream" 
that at an earlier time was described as "downstream" . In this sense, the pre- 
sented nonlinear ray equations are valid for weak shock waves, if there are no 
subsequent shock front intersections. Nevertheless, to investigate plasma heating 
we calculate the caustic pattern, neglecting the modified downstream plasma, 
which is correct for weak shock waves. We note again that since subsequent 
evolution of the downstream plasma flow cannot be analysed with the method 
used, we do not discuss the transformation of a null point into a current sheet 
under the shock wave effect. 

For a linear wave approaching the null point in a cold plasma, the analytical 
law of the amplitude growth has been found. The wave amplitude increases 
exponentially and therefore the initially linear wave transforms into a shock, 
which has encouraged us to investigate the case of shock wave propagation. 
This analytical law expressed by Equation (and Equation |QJ) is consistent 
with the analytical solution that was obtained for the problem of the incidence 
of a cylindrically symmetric fast-mode wave on a 2D null point ( |Syrovat skii, 
ITMrJl |Craig and McClymontl ITMT1 |Longcope and Priest) 12007)) . Indeed, in their 
solution the radial component of the plasma velocity in the wave grows as 1/r 
with the distance r to the null point. In the cylindrically symmetric case all the 
rays are straight and directed to the null point. So we have the following ray 
equation: dr/dt = —Va = —r, and hence we obtain: r oc e , i.e., the wave 
amplitude grows exponentially, like that in Equation ([3]). The additional time 
dependence in Equation ((3|) is due to the plain (not cylindrical) initial wave 
front and represents the fact that all the rays are different and depend on initial 
coordinate xq. 

A complex caustic is formed around the magnetic null point. The wave ampli- 
tude tends to infinity at the caustic, resulting in plasma heating. The heating is 
distributed in space and occurs mainly out of the Va = c layer if the initial wave 
front is rather extensive. However, the first heating associated with a central part 
of the wave occurs within the Va = c layer. On the contrary, in the well-known 
case of a linear wave in a cold plasma, plasma heating takes place in a small 
vicinity of the null point. Since real null points in the solar corona appear to 
capture only a (central) part of the wave, while its outer part propagates away, 
the size of the Va = c layer may be chosen as the spatial scale of plasma heating. 

Nonlincarity of a wave is crucial for its propagation. A shock wave is able 
to pass through a null point even in a cold plasma. The growth of the shock 
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wave amplitude is restricted, since shock waves decay significantly due to the 
nonlinear energy dissipation in the shock front. Therefore along with the heating 
at a caustic, a substantial heating of plasma occurs in the neighbourhood of the 
null point. A significant part of the shock wave energy transforms into heat there. 

Drawing the caustic pattern, we have explained the complicated surface of 
the wave front. A break of the wave front occurs at the caustic, thus triangles of 
the wave front arise when the front passes through a caustic cusp. The vertices 
of the triangles travel along the caustic lines. 

To summarise we discuss a chance to estimate the magnitude of plasma heat- 
ing associated with propagation of waves near a magnetic null point. Conversion 
of the wave energy into heat occurs due to dissipative processes in a plasma when 
the gradients of MHD quantities increase significantly. In a linear wave the gra- 
dients grow owing to the wave extent decrease as well as the amplitude growth as 
the wave approaches a null point. To estimate the value of the energy converted 
into heat in this case, one should consider non-ideal magnctohydrodynamics. 
In turn, in a shock wave the gradients of MHD quantities also originate at the 
wave front due to the nonlinear steepening of the wave profile. Here, the value of 
the energy converted is taken into account by the laws of shock wave damping. 
Since the sound speed in the real solar coronal plasma does not drop to zero, 
the wave extent can remain rather large, and so viscosity, thermal conductivity 
and finite plasma conductivity do not result in energy conversion. In such a 
case considerable plasma heating will occur due to the wave energy dissipation 
in the shock front, with the nonlinear heating predominating over heating at 
caustics. So, evaluating the difference between the shock wave amplitude and the 
amplitude calculated without dissipation in the discontinuity (as in Figure [6]), 
one can provide some rough estimates about the plasma heating magnitude. 



Appendix 

A. Ray equations 

The sets of ray equations have the following form in Cartesian coordinates: 

dx _ ( N KU s h \ k x ^ da dk x _ da ^ 

dt \ 2 J k dk x ' dt dx ' 

dt [ a + N 2 J k +k dk z > dt ~ dz K (5) 

where k x , k z and k are the components of the wave front normal k and its 
magnitude respectively, and N the numerical factor switching the linear and 
nonlinear regimes of the wave propagation. We should take N = for studying 
a linear wave and N = 1 for a weak shock one. 

The adjoint set consists of differential equations for the derivatives dx/dr], 
dz/drj, dk x /dr] 7 dk z jdr\ and is derived by differentiating Equations ([5]) (with 
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N = 0) with respect to the ray coordinate r\: 

d ( dx \ k x da a dk x ak x dk da dk d ( da 

dt \ dij J k drj k drj k 2 d-q dk x drj dr\ \ dk x 

d / dz\ k z da a dk z ak z dk da dk ^ d ( da 

dt \ dr] J k drj k drj k 2 drj dk z drj dr\ \ dk z 

d ( dk x \ dk da ^ d / <9a\ 



dt \ drj J d-q dx drj \ dx 
d ( dk z \ dkda ^ d /da 
dt \ drj J drj dz d-q \ dz 



where 



da s-^ ( da dr a da dk a \ dk 1 dk a 
dr] ^ \dr a drj dk a drj J drj k ^ a drj 

d ( da \ >— ^ / d ( da \ dr a d ( da \ dk a 
dr] \dkp J " \dr a \dkp J drj dk a \dkp J dr] 

d ( da\ { d ( da\ dr a d ( da\ dk a 

dr) \drp) ^ \dr a \drp) d-q dk a \drp) d-q 
r a x, z, k a — k Xj k z . 

The numerical coefficient k depends on the plasma beta as well as the angle 
a between the wave front normal and the magnetic field as 

(1 + p + Qi)(/3( 7 - 2) + 3Qi) - 2^ cos 2 q( 7 - 1 - + Qi) 
2Qi(l + /3 + Qi) 



where Qi = + fi) 2 — 4/3 cos 2 a, j3 = c 2 /Vj, and 7 = 5/3 is the adiabatic 
index. Values of k are restricted by the limits 1/2 < k < 3/2 and are shown 
in Figure |S] When the wave crosses the Va = c layer along the magnetic field 
(a = 0), k has a jump. 

Set of ray equations ([5]) is not closed in the nonlinear case since it includes the 
wave amplitude U s h- It can be calculated with the laws of solitary shock wave 
damping (jUralovl I1982p , which are derived by using values of the amplitude and 
duration of a simple wave from which the shock forms: 

Ti \ dT\ KUl 



where T\ is the duration increment of the simple wave with amplitude u\\ Tq is 
the initial duration of the shock wave. The value of u\ can be found from an 
equation similar to Equation @. 

Solving sets of ordinary differential equations ([5]) , (O and ([7]) numerically we 
are able to model the propagation of fast-mode MHD waves. 
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Figure 8. The values of coefficient ft. The points at the surface denote the relevant values 
corresponding to some particular cases of propagation. There is a jump of re at a = and 



B. Derivation of nonlinear ray equations 

The sets of ray equations can be obtained with a different method which does not 
use the eikonal expansion of the MHD equations. In fact, this method is different 
only in form and it is more convenient for analysing nonharmonic waves which 
we deal with in this paper. We apply this method for deriving the set of nonlinear 
ray equations. 

The moving surface of the wave front can be described by the equation 
^(r,^) = 0, or in the differential form d<f> = 0. Let a small surface element 
displace by dr in time dt. Then we have: 



where dr/dt = q is the velocity of the front element and grad$ the wave front 
normal by definition. The scalar product can be rewritten as q n |grad$|, with 
the q n denoting the normal speed of the wave front. So we obtain a partial 
differential equation. 

If we consider a linear fast-mode MHD wave, q n = a, and we have the same 
"eikonal" equation and the same set of ray equations as presented in Section EOl 
and Appendix [Al (the case of N = 0), respectively However, we are interested 
in the weak shock wave propagation. The normal speed of the weak shock front 
can be derived from the jump conditions across a MHD discontinuity, and it is 



13=1. 




(8) 
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equal to a+K,U s h/2 in the nonlinear acoustics approximation. Substituting it into 
Equation ([5} and solving with the method of characteristics, we obtain the set of 
ray equations, which is similar to the linear one and in which substitution a — > 
a + KU s h/2 has been performed. In these new ray equations the additional terms 
dU s h/dk a and dU s h/dr a appear. We neglect them, keeping only quantity U s h 
itself and omitting its derivatives, exploiting the nonlinear geometrical acoustics 
approximation. Thus we obtain the set of nonlinear ray equations ([5]) (with 
N = 1) presented in Appendix |A"1 
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